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TRIDIAGONALIZATION OF THE HYPERGEOMETRIC 
OPERATOR AND THE RACAH WILSON ALGEBRA 

VINCENT X. GENEST, MOURAD E. H. ISMAIL, LUC VINET, AND ALEXEI ZHEDANOV 


Abstract. The algebraic underpinning of the tridiagonalization procedure is 
investigated. The focus is put on the tridiagonalization of the hypergeomet¬ 
ric operator and its associated quadratic Jacobi algebra. It is shown that 
under tridiagonalization, the quadratic Jacobi algebra becomes the quadratic 
Racah-Wilson algebra associated to the generic Racah/Wilson polynomials. 
A degenerate case leading to the Hahn algebra is also discussed. 


1. Introduction 

This paper aims to offer an algebraic framework for the tridiagonalization ap¬ 
proach to the study of orthogonal polynomials. This procedure has been inves¬ 
tigated by Ismail and Koelink in a series of papers [101 HU IE]- In [El Section 
3], these authors have shown in particular that a special case of the Wilson poly¬ 
nomials can be obtained from the Jacobi polynomials by tridiagonalization. The 
hypergeometric polynomials of the Askey tableau can be described from the most 
general to the more specific by starting with the 4-parameter Wilson polynomials, 
which sit atop the tableau m Chapter 9]. It is however appealing to develop an 
understanding of the Askey scheme from the bottom up by constructing the more 
involved families of polynomials from the simpler ones. This is in part what the 
tridiagonalization method allows to do. The idea of “lifting” within the Askey 
scheme is also in that spirit Sections 15.1, 15.2]. 

The main concept behind the tridiagonalization approach is to start from an 
operator L that has a family of orthogonal polynomials {Pn{x)} as eigenfunctions 
and to combine it with other operators so as to obtain a new operator M that 
acts on {Pn{x)} in a tridiagonal fashion. The diagonalization of the operator M 
by another family of orthogonal polynomials {( 5 n( 2 ;)}, when possible, will relate 
{P„(a;)} and {Qnix)}, whose characterization will benefit from the setup. 

It is known that many of the important properties of the classical orthogonal 
polynomials are encoded in polynomial algebras and their representations [15] . This 
is done for instance by using as generators the two operators involved in the bispec- 
trality of the family of polynomials of interest. We here integrate the tridiagonaliza¬ 
tion method into this algebraic picture. We focus on the connection between Jacobi 
and Wilson polynomials via tridiagonalization and explain how the Wilson/Racah 
algebra emerges in the process. 

The outline of the paper is as follows. The tridiagonalization of the hyperge¬ 
ometric operator L is revisited in Section 2. The approach of m is generalized 
by considering both left and right multiplication by x. It is observed that L and 
the tridiagonalized operator M interchange their role under the transformation 
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X x~^ . In Section 3, it is shown that the generic Wilson polynomials with 4 
parameters arise from the tridiagonalization of L. This leads naturally to Koorn- 
winder’s integral representation for the Wilson polynomials M- The case when 
both L and M preserve the space of polynomials of degree smaller or equal to a 
positive integer N is considered. This finite-dimensional reduction is the object of 
Section 4. The connection with the Wilson/Racah algebra is explained in Section 
5, where it is shown that the hypergeometric generator L and the tridiagonal- 
ized operator M realize this algebra. Reference is made to the realization of the 
Wilson/Racah algebra out of the intermediate Casimir operators occurring in the 
coupling of three irreducible representation of su(l,l). The operators introduced 
in the tridiagonalization of L will be observed to coincide with the Casimir opera¬ 
tors arising in the one-dimensional model of the Racah problem elaborated in [6]. 
Section 6 is dedicated to the degenerate case where the parameters of the tridiag- 
onalized operator M are such that it becomes a first order differential operator. It 
is then seen that the Hahn algebra and polynomials occur. 


2. Tridiagonalization procedure for the hypergeometric operator 


In this section, the tridiagonalization procedure for the hypergeometric opera¬ 
tor is revisited and generalized. The quadratic Jacobi algebra associated to this 
operator is also presented. 


2.1. The classical Jacobi polynomials and the quadratic Jacobi algebra. 

Consider the ordinary hypergeometric operator 

(2.1) L = a;(l — x)d‘l -\- (^a -\-1 — [a + (3 + 2) x^dx, 

where a,/3 are real parameters satisfying a, 13 > —1. The operator L possesses an 
infinite family of polynomial eigenfunctions that obey the eigenvalue 

equation 

(2.2) LP^«(x)=A„Pi“>«(x), 
with eigenvalues 

(2.3) Xn =-n{n + a + P + 1). 


The monic polynomials 


(2.4) (x) = ( + f 

(a-|-p-|-n-|-l)n \ 


—n n + a + (3 + \ 
q; -I- 1 


; X 


with {a)n denoting the Pochhammer symbol 

(a)„ = a{a + l)(a -I- 2) • ■ • (a -I- n — 1), (a)o = 1, 
coincide, up to an affine transformation, with the Jacobi polynomials [T]. The 
polynomials (x) satisfy the three-term recurrence relation 

(2.5) Pi'^+Pix) + 6„Pi“’«(a;) + u^P!rJ\x) = xP^^^^^x), 

with the recurrence coefficients 

n(n + a){n + I3){n + a + j3) 


— 


— 


(2n + a + /3 — l)(2n + a + /3)^(2n + a + /3 + 1 ) ’ 

1 ( 1 _ 1 

2/1 -\- o. -\- P ‘2 ti -\- ol -\- P -\- ‘I 
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They obey the orthogonality relation 


( 2 . 6 ) 


X°‘{1 — x)^dx = hn Sn 


where the normalization factor is given by 


(2.7) 


r(a + i)r(/3 + i) 
= r(o + /^ + 2 ) 


where r(a;) is the Gamma function. 

There is an associative algebra known as the quadratic Jacobi algebra which 
is associated to the hypergeometric operator and to the Jacobi polynomials [SI 
Formula (3.4) and Section 4]. This algebra can be obtained as follows. Let Ki and 
K 2 be defined as 


(2.8) iFi = L, K2= X, 

where L is given by (EH) and where X is the “multiplication by a:” operator; that 
is X f{x) = xf{x). Introduce a third operator K 3 , defined as K 3 = [Ki,K 2 ], where 
[a, &] = ab — ba is the commutator. From (12.11) . it follows that K 3 has the expression 

K 3 = —2x{x — l)dx + (l + a— (a + /3 + 2)a;). 

By a direct calculation, one finds that iFi, K 2 and K 3 satisfy the relations 

[K3,K2]=K3, 

(2.9) [K2,K3]=a2Kl + dK 2 , 

[K 3 , Ki] = a2{Xl, K 2 } + diFl + 02^2 + 62 , 

where {a, b} = ab ba and where the structures constants 02 , d, C 2 , 62 are given by 
02 = 2 , d = — 2 , C 2 = —{a + /3)(q! + /3 + 2 ), 62 = (o + 1 )(q! + /3). 

The relations (12.91) define the Jacobi algebra. This algebra has a Casimir operator 

(2.10) Q = a 2 {Kl Ki} + Kj + (a^ + c 2 )iF| + 47^1,7^2} + {da 2 + 2e2)K2, 
which commutes with all generators Ki, K 2 , K 3 . In the realization (|2.8I) . one has 

(2.11) Qf{x) = {a^ -l)f{x). 

The polynomials form a basis for infinite-dimensional irreducible 

representations of the quadratic Jacobi algebra (12.91) . In this basis, the operator 
Ki is diagonal with eigenvalues EH, and the operator K 2 is tridiagonal, as is 
manifest from ( 1 ^ . The monomials also form a basis for these irreducible 

representations. In this basis, both generators Ki and K 2 are two-diagonal. Indeed, 
one has 


(2.12) iFi x” = L x" = A„a;" + with €n = n{n + a). 


and evidently X x^ = 
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2.2. Tridiagonalization of L. Consider the operator 
(2.13) M = tiXL + T 2 LX+ T 3 X+ To, 

where r^, i = 0,1, 2, 3 are arbitrary parameters such that T 1 +T 2 7 ^ 0. Without loss 
of generality one can take n + r 2 = 1 , as this corresponds to scaling the operator 
M. Upon using (EH) in (12.131) . one finds that M has the expression 


(2.14) M = x^{l — x) d1 + X {a + I+ 2t2 — {a + fi + 2 + 2t2)x) dx 

- {T 2 {a + /3 + 2) - T 3 ) a; + (a + 1 )t 2 + tq. 


Remark. The condition ri + T 2 yf 0 is essential for the operator M to be a second- 
order differential operator. If ti -|- r 2 = 0 , M becomes a first-order operator. 


From (12.211 . (12.511 . (12.1211 and (12.1311 . it is seen that the operator M has the 
following properties: 

• It acts in a three-diagonal fashion on the polynomials Pn°‘’^\x). One has 
indeed 


(2.15) MPi^’P\x) = (riA„ + T 2 A„+i + T^)Pi"^^\x) 

+ i^nbn + T3bn + tq)P^’^\x) + (riA„ -|- T 2 A„_i -|- T 3 ) P^_’f^(a:). 

• It acts in a two-diagonal fashion on the monomials x”. One has 

(2.16) Mx" = (riA„ -f T 2 A„+i -b T 3 ) x"+^ -b (rie„ -b r 2 e„+i -b tq) x", 
where A„ and e„ are respectively given by (12.311 and (12.121) . 

There is an interesting duality between the operators L and M. Indeed, one ob¬ 
serves using the explicit expressions ( 1 ^ and (12.131) with the condition n -b T 2 = 1 
that the operator L can be expressed in terms of the operator M as follows: 

(2.17) L = TiX-^M + T2MX-^ - (to -b 2tiT2)X-^ -b 2 tiT2 - T 3 , 


where X~^ is the “multiplication by x“^” operator, i.e. X~^ /(x) = x~^ f{x). The 
relations (12.131) and (12.1711 have the same structure under the substitution X 
X~^. Consequently, one can expect that the operator L will be three-diagonal in 
an appropriate (discrete) basis of eigenfunctions of M. It is convenient to introduce 
two new parameters vi and 1/2 instead of ri and T 3 


(2.18) T 2 =-( 1 -b 1^1 - 1 ^ 2 ), 


T3 = + 


(a -b/3)(i^i -b 1^2 - 1) 


Consider the change of variable x = Ijy together with the similarity transformation 


(2.19) L = y''^-^ M = y’'^-^ M 


Under these transformations, the operator L takes the form 

(2.20) L = y'^iy - l)dy + y{aiy + bi)dy + Ciy -b dx, 
where the coefficients oi, 61 , ci, dx are given by 

ax = 3 — a — 21 ^ 2 , bx = a + f3 + 2z/2 — 2, 

Cl = (z^2 - 1 )(q; -b 1^2 - 1 ), dx = —{1^2 — 1 )(q: + 13 + 1^2)- 
The operator M becomes 

( 2 . 21 ) M = y{y - l)dl + {a 2 y + b 2 )dy + d 2 , 

















TRIDIAGONALIZATION AND THE RACAH-WILSON ALGEBRA 


5 


where a 2 ,& 2 ,(i 2 read 

02 = 2 — — ^2 — a, &2 = a + /3 + — 1, 

d2 = To + viV2 + - l){yi + V2- !)■ 

As is seen from (|2.20l) and (12.211) . the operators L and M have their roles inter¬ 
changed under the change of variable y = Ijx and the similarity transformation 
(I2.19|) . Indeed the operator M coincides with the hypergeometric operator (up to 
an affine transformation), while L has a form similar to M. This means in partic¬ 
ular that there exists a set of eigenfunctions {V'n(a^)} for n = 0,1,2,... that satisfy 
the eigenvalue equation 

M 1pn{x) = Xn-ipnix), 

and on which the operator L acts in a three-diagonal fashion 

( 2 . 22 ) Llpnix) = P^n'lpn+l{x) + Tjn'lpriix) + C,npJn-l{x), 

for some coefficients ^n,^ra,Cn- Indeed, since M is the ordinary hypergeometric 
operator, it possesses a set of polynomials as eigenfunctions. It is verified that M 
admits the discrete set of eigenfunctions 

(2.23) ipn^x) = x'''^~^P^°‘’^\l/x), Xn = n{n - 1) + a 2 n + d 2 ^ 

where /3 = /3, a = —a — jd — vi — V 2 - The property (12.221) then follows from the 
fact that Pn°‘’^\x) are the polynomials (12.41) satisfying the three-term recurrence 
relation (12.51) . 

3. Wilson polynomials arising from tridiagonaliation 

In this section, we show that the tridiagonalization of the hypergeometric oper¬ 
ator L leading to the operator M gives rise to the generic Wilson polynomials. 

Let ip{x\ A) be an eigenfunction of M corresponding to the eigenvalue A 

(3.1) Mil){x] A) = Ktpix] h), 

which is integrable with respect to the Jacobi weight. Consider the expansion of 
^(x; A) in the formal Jacobi series 

OO 

(3.2) i/>(x;A)=^Gfe(A)Pi“’«(x). 

Let us write the expansion coefficients as 

Gfe(A) = Go(A)Sfe Qfc(A), 

where S*, are normalization coefficients, (5o(A) = 1 and where Go (A) is given by 

I 

Go(A) = — / 'il){x] k)x^{l — xYdx. 

do Jo 

Upon applying M on both sides of (13.21) and using (12.151) . one finds that the coef¬ 
ficients Qk{A) satisfy the three-term recurrence relation 

Aii)Q„+i(A) + eY^QYA) + eY'>Qu-i{A) = AQ„(A), 












6 


V.X. GENEST, M. E. H. ISMAIL. L. VINET, AND A. ZHEDANOV 


where 

EI^'> = + T2Xn + Ts), E^'> = X^bn + T^hn + Tq, 

‘n 

E^) ^ ::^zl(^riXn-l + T 2 Xn + T 3 ). 

“n 

Let us choose the coefficients 5„ to ensure that 

E^-^ = + T2Xn + T 3 ) = 1, 

and that Sq = 1. These requirements determine S„ uniquely. One then finds that 
Qn{A) are polynomials in A with the three-term recurrence relation 

(3.3) gn+l(A) + Bn Qn{A) + Un Qn-l{A) = A Q„(A), 
recurrence coefficients 

Bn=E^'> = Xnbn + Tsbn+To, C/„ = (ti A„_i-|-T 2 A„-|-T 3 ) (ti A„-|-r 2 A„_i-|-Ta) , 

and initial conditions (3_i(A) = 0 and Qo(A) = 1. It follows that Q„(A) are 
orthogonal polynomials in A. These polynomials can be identified with the Wilson 
polynomials W„(a:; oi, 02 , 03 , 04 ), which depend on 4 parameters 04 , 01 , 03,04 and 
which are defined by the three-term relation 

(3.4) W'„+i(a;) + (A„ + C7„ - aj) W„(a;) + A„_iC„W„_i(ai) = xW'„(ai), 
where 

^ {n + g — l)(n -I- oi -I- a 2 )(n -I- oi -I- a 3 )(n + ai + 04 ) 

(2n + 5 -l)( 2 n + g) ’ 

_ n(n + 02 + as — l)(n + as + a 4 — l)(n + as + 04 — 1) 

" (2n-I-(? — 2)(2n-I-g — 1) ’ 

and where g = oi -I- 02 -I- 03 -I- 04 . Upon comparing ()3.3I) and (13.4L one finds 

(3.5) Q„(A) = IT„( 7 -A;oi,02,a3,a4), 
with 7=1 (oi 4-02 — 04 — 02 - 1 - 2 to) and where 

a = oi -I- 02 — 1, /3 = 03 -I- 04 — 1, 

1^1 = 1 — Oi — O 3 , 1^2 = 1 — 02 — 03 . 

We have thus identified the polynomials Qn(A) with the generic Wilson polyno¬ 
mials. This indicates that the tridiagonalization procedure, when applied to the 
hypergeometric operator L, leads to the generic Wilson polynomials, which depend 
on four parameters. 


Remark. While four parameters appear in the definition (12.131) of M, only two of 
them are essential since tq is merely a translation and since M can be scaled by 
an arbitrary factor. Together with the two parameters a, /3 of the hypergeometric 
operator, this accounts for the four arbitrary parameters of the Wilson polynomials 
arising in (13.51) . 


Let us parametrize the eigenvalues A appearing in (13.11) 

. 1 — a, ,, , (i^i - r'2)(Q;-I-2g-I-1) 

A = To -f -f (a -b q){q -b 1) -b ^^1 - L, 


as 


2 


2 
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where g is a new parameter. It is verified that dSU admits the solution 
(3.6) V'(a;;A) = x)~^ 

where ai, /3i and 71 are given by 

ai = a + I + vi + q, j3i = I — (3 — V 2 + q, "fi = 2 + 2q + a + vi — V 2 . 

Note that one has 71 — oi — = /3 and hence it follows that for g > 0 and /3 > 0 

the integrals 

f x'^{l — x)^x'^tf}{x;A)dx, 

Jo 

exist for all non-negative integer n. Thus the function ip^x^h) taken as in (13.61) 
satisfies the required property of being integrable with respect to the Jacobi weight 
function. The other independent solution of (j3.1ll does not have this property and is 
hence discarded. Formally using ()3.2I1 suggests the following integral representation 
of the generic Wilson polynomials (5„(A): 

(3„(A) = —— f 'ii){x]h)P^°''^\x)x^{l-xYdx. 

nn'~'0V^)^n Jo 

This formula is seen to coincide with Koornwinder’s integral representation for the 
Wilson polynomials [M]. It is interesting to note that the dual solutions tpn{x) 
described by (12.231) do not satisfy the condition that the integrals 

f ljjnix)x°‘{l — x)^dx 
Jo 

exist for all non negative integers n. Indeed, the function ipnix) has the polynomial 
part Pn°'’^\l/x). Hence, for sufficiently large n, it is seen that the integral 


[ x-^x°^{l-xfdx 

Jo 


will diverge. 


Remark. In the process of describing the algebraic underpinning of the tridiago- 
nalization procedure, we have seen the integral representation of the Wilson poly¬ 
nomials of m arise as a side result. As pointed out by a referee, one could have 
started instead from that integral representation to derive the embedding of the 
Racah-Wilson algebra into the Jacobi algebra. 


4. Finite-dimensional reduction 

In general, the two discrete sets of eigenfunctions {Pn°'’^\x)} and {■)/;„(a:)} of 
the operators L and M are infinite-dimensional. However, there is a case for which 
it is possible to work with finite-dimensional sets. This is the object of this section. 

Consider the action of the operators L and M on the monomials x". In view of 
(I2.I2L it is manifest that L preserves the space of polynomials of a given degree 
N, where iV is a non-negative integer. In contrast, it is visible from p.l6|) that in 
general M does not preserve the space of polynomials of a given degree, as it maps 
polynomials of degree n to polynomials of degree n -I- 1. Consider the special case 
for which the condition 

( 4 . 1 ) tiAat -I- T2AJV+1 + 7-3 = ri(A 7 v — Aat-i-i) + Avr+i -I- T3 = 0 , 
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holds for some non-negative integer N. It is clear from (12.161) that in this case the 
operator M preserves the space of polynomials of degree less or equal to N. Thus 
under condition 63) both operators L and M preserve the space of polynomials 
of degree less or equal to N and one can restrict the action of these operators to 
the N + 1-dimensional space spanned by these polynomials. There are two natural 
bases for this space. The first one has the polynomials Pn°‘’^\x) for n = 0,1,..., 
given in (1^ as basis vectors; this is the eigenbasis of L. The second one has 
the functions ipn{x) ioi n = 0,1,..., N given by (12.231) as basis vectors; this is the 
eigenbasis of M. In the parametrization (I2.18|) . the truncation condition (14.11) reads 

{v 2 -N - 1){N +l + i^i+a + P) =0. 

It follows that one can take either 1/2 = A^-l-l ori^i = —N — 1 — a — p. For 
definiteness, we choose to use the condition 1/2 = N + 1. This leads to the following 
expression for the basis elements ipn{x)\ 

(4.2) V'n(cr) = x^ 

It is easily seen that these are polynomials of degree N — n. Consider the expansion 
of the basis functions tp„ {x) given by in terms of the basis functions 
one writes 

N 

^n{x) ='^RnkPt’^\x). 

k=0 

The expansion coefficients Rnk can be obtained from the orthogonality relation 
(12.61) of the basis functions P^'^\x). One finds 

Rnkhk = / X"(l - xYdx, 

Jo 

which can be considered as a finite-dimensional analog of the Fourier-Jacobi trans¬ 
form proposed by Koornwinder |14) . In view of (13.51) . this formula leads to an 
integral representation for the truncated Wilson (Racah) polynomials that arise in 
the expansion coefficients Rnk- 

5. TrIDIAGONALIZATION AND THE RACAH-WiLSON ALGEBRA 

In this section, the relation between the tridiagonalization of the hypergeometric 
operator, the Jacobi and Racah-Wilson algebras is discussed. It is shown that 
the operators L and M provide a realization of the Racah-Wilson algebra. The 
connection between this realization and the one arising from the Racah problem 
for the Lie algebra su(l, 1) is also established. 

5.1. The Racah-Wilson algebra. The Racah-Wilson algebra [5] is generated 
by the elements Ai, A 2 and their commutator A 3 . These elements satisfy the 
commutation relations 

[^ 1 ,^ 2 ] = A 3 , 

[A2, A3] = ai{Ai, A2} -l- 02^2 -|- 7 iAi -|- 5A2 -b ei, 

[A3, Ai] = a2{Ai, A2} + aiAi -j- 72A2 -I- SAi -b £2, 


(5.1) 
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where ai, a 2 , 71 , 72 , ei, £2 are real structure constants. The algebra (lEID has a 
Casimir operator which has the expression 

R = ai{Af, A 2 ] + Q!2{^2i ^ 1 } + (o^i + 7i)^i + (q^2 + 72)^2 + ^3 

+ (5 + aia 2 ){Ai, A 2 } + {Sai + 2 ei)Ai + {8012 + 2 e 2 )^ 2 ) 

and which commutes with all generators Ai, A 2 , A 3 . Upon comparing the defin¬ 
ing relations dsn of the Racah-Wilson algebra and those of the quadratic Jacobi 
algebra (j2.9ll . it is directly observed that the latter is a special case of the former 
that corresponds to taking ai = = ex = 0. The Racah-Wilson algebra is the 

algebraic structure associated to the Racah and Wilson polynomials, as it can be 
realized in terms of the operators involved in the bispectrality of these families of 
orthogonal polynomials. This realization, which can be found in for example, 
is similar in spirit to the realization of the Jacobi algebra (12.911 in terms of the 
operators L and X exhibited in Section 2. For a review of this algebra and of some 
of its applications, one can consult [?]■ 

5.2. Racah Wilson algebra arising from tridiagonalization. We now show 
how the Racah-Wilson algebra arises in the present context. Let Ai{x) and A 2 {x) 
be defined as 

(5.2) Ai{x)=L, A 2 {x) = M, 

where L and M are as in (12.111 and (12.1311 . respectively. Moreover, let ^ 3 ( 2 ;) read 

A^ix) = [Ai{x),A 2 {x)\. 

Then, by a direct calculation, it is verified that the operators ^ 1 ( 2 :), ^ 2 ( 2 ;), ^ 3 ( 2 ;) 
satisfy the defining relations dSU of the Racah-Wilson algebra with the structure 
constants taking the values 

oi =-2, 02 = 2, (5 = a(a-I-/3-I-1)-f/3 - 4ro - 2 t3, 

71 = 1 - 4(ro +T 2 - r|), 72 = -(o -h /3)(a + 15 + 2), 

ei = 2 tq +{1- a^)T 3 -h 2 x 3 X 3 2(a l)(/3 -f l)(x| - X 2 ) - {a + l)(a -h /3)xo, 

€2 = (a + /3)((a + (5 + 2)xo + (a -I- l)x 3 ). 

In this realization, the Casimir operator acts as a multiple of the identity 

R{x) = 2{a + l)(/3 -h l)(a /3)(xf - xi) + {{a + P + 1)^ - 5)xq 

-f 2(a -h l)(a -f /3 )(xoX 3 -h xq + X3) - 4x3X3 -h (o^ - l)x|. 

In terms of the generators (12.811 of the quadratic Jacobi algebra ( 1 ^ . the generators 
(I5.2[l of the Racah-Wilson algebra have the expression 

^1(2:) = Ki, A 2 {x) = xiRr 2 Ki -I- T2K1 K2 + T3K2 + xq. 

This provides an embedding of the Racah algebra into the Jacobi algebra in their 
standard realization on the space of polynomials. 

Remark. In lieu of M, one can introduce an operator J by applying the tridiago¬ 
nalization procedure to the hypergeometric operator L using instead the “multipli¬ 
cation by 2 ; — 1 ” operator; this operator is defined by means of 

(5.3) J = p^{X- l)L + P2L{X - 1) + p3{X - 1) + p 3 , 
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where pi, p 2 , ps, po are arbitrary parameters such that pi + p 2 ^ 0. Due to the 
symmetry of the hypergeometric operator with respect to the change of variable 
a; —>■ 1 — a;, it is clear that the operator J has the same properties as the operator 
M. In particular, one obtains the generic Wilson polynomials when expanding the 
eigenfunctions of the operator J in a series of Jacobi polynomials. Furthermore, 
it is verified that if one defines the operators Bi{x) = L, and B 2 {x) = J as well 
as B^{x) = i? 2 (a;)], one obtains another realization of the Racah-Wilson 

algebra (ED- 


5.3. Connection with the Racah problem. There is an interesting interpreta¬ 
tion of these results in terms of the Racah problem for the su(l, 1) algebra. Consider 
the su(l, 1) algebra with generators Sq, S+, S- and commutation relations 

[So,S±] = ±S±, [S-,S+]=2So. 


The Casimir operator 

C = S^-So- S+S-, 


commutes with all generators. On the irreducible representations of su(l, 1) that 
are relevant to the present paper, the Casimir operator takes the value cr(cr — 1). 
For cr > 0, the irreducible representations of the positive-discrete series are defined 
on the space spanned by the basis vectors e„, n = 0 , 1 ,..., by the following actions: 


SqG-ji — (tZ J- (7^€-yi, — 'y^j^iCn+l, S—€n — 'Yn^n — l: 


where 7 „ = ^Jn{2a -I- n — 1). Consider the tensor product of three representations 
of the positive-discrete series with representation parameters cti, CT 2 and 173 . We 
use the notation 



c(i) 

*^0 






hj = 1,2,3, 


where the superscript i in 5'^*^ specifies on which representation space the generator 
acts. The intermediate Casimir operators Cik are defined as 

ak = for (ifc) e {(12), (23), (13)}, 

and the total Casimir operator C is written as 

Q = _ 5 - 1123 ) _ ^(123)^(123) 

where 5^1 = 5^1 -|- + 5'^ ^ and similarly for . By construction, the total 

Casimir operator commutes with all intermediate Casimir operators. 

Consider the decomposition of the three-fold tensor product representation in 
irreducible components. On each of these components, the total Casimir operator 
takes the value CT 4 (cr 4 — 1), where 174 = + tri + CT 2 J- us and is a non-negative 

integer. For a given N, there corresponds A^ J- 1 eigenvalues of the intermediate 
Casimir operators Cy which are of the form crij{aij — 1) where aij = kij + ai + aj 

}12') f23^ 

and where kij = 0 , 1 ,..., N. Let and be eigenbases for C12 and C23, i.e. 

Ci2e]^^^ = ai2iai2 - l)e[,^j\ = ^ 23 ( 1^23 - 1)6^^, 

where fci 2 , ^23 = 0,1,..., A^. The Racah problem consists in obtaining the expan¬ 
sion coefficients of the basis over the basis 

K23 ^12 
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For more details on the Racah problem for su(l, 1), the reader is referred to [4]. 

A one-dimensional model for the Racah problem was constructed in [ 6 ] by sep¬ 
aration of variables and dimensional reduction. It allows to express the three in¬ 
termediate Casimir operators as second-order differential operators in one variable. 
One has 

012 = - x)d'^+ x[{N - 1 - 2ai)x + 2{ai + a 2 )]dx 

+ 2NaiX + (di -I- (T2)(cri -I- tT2 — 1 ), 

4 023 = x(x -l)d^ + [2(1 - N - a 2 - a 3 )x -I- (iV - 1 -|- 2 a3)\dx 

-I- (Af -I- (73 -I- a 2 ){N + (73+ a 2 - 1 ), 

O 31 = x{x — 1)^31 -b (1 - x)[{N - 1 - 2ai)x + 1 - N - 2 a 3 ]dx 
+ 2Nvi{\ — a;) -b ( 0-3 -b cri)(cr 3 -b di — 1). 

Any pair of operators Cik satisfy the Racah-Wilson algebra [ 6 ] . Up to affine trans¬ 
formations, the operators L, M and J introduced in (12.1|) . (12.1311 and (15.3|) . coincide 
with the operators Cik of the one-dimensional model (15.4|) . This provides an expla¬ 
nation for the fact that these operators satisfy the Racah-Wilson algebra. More¬ 
over, it shows that the tridiagonalization procedure of the hypergeometric operator 
is intimately related to the Racah problem for su(l, 1). 

6 . Degenerate case. The Hahn algebra 

In this section, we consider the degenerate case for which ri -b 12 = 0 in the 
tridiagonalization procedure. It is shown that in this case the Hahn algebra and 
Hahn polynomials arise. 

Let Ti -bT 2 = 0 in (12.131) . One can set n = 1 and T 2 = —1. With these conditions, 
the operator M has the expression 

(6.1) M = 2x{x — l)dx -b (a -b /3 -b T 3 -b 2)a; -b tq — a — 1. 

It is observed that M is now a first-order differential operator. Let 

(6.2) T 3 = -2N -2-a- P, 

with N a non-negative integer. It is seen that both operators L and M preserve the 
N + 1-dimensional space of polynomials of degree less or equal to N. It is already 
known that the eigenfunctions of L are the polynomials ()2.4I) . The polynomial 
eigenfunctions ipnix) of the operator M have the expression 

(6.3) V'n(x)=cc"(l-a;)'^-", 
and satisfy the eigenvalue equation 

(6.4) MV'n(x) = (tq - 2n - a - l)^/>„(a:), 

for n = 0,1,..., A^. By construction, the operator M is 3-diagonal on the polyno¬ 
mial basis PiC‘'^\x)\ 

(6.5) MPt^\x) = (A„ - A„+i + T 3 )P!tpf\x) 

+ (t 3&„ -b To)P^°‘’^\x) + u„(A„ - A „_1 -b t 3 )P^°P^\x), 
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where A„ is given by (j2.3l) . It is easy to show that the operator L acts in a tridiagonal 
fashion on ipnix). Indeed, one has 

(6.6) L-ijjnix) = {N - n){N - n +/3)^n+iix) 

+ (2n^ + (a — /3 — 2N)n — {a + l)N)'ipn{x) + n{n + 

Applying the same method as in Section 2, one arrives at orthogonal polynomials 
(5n(A) satisfying the three-term recurrence relation (13.31) with coefficients 

(6.7) Bji — To, Un — Xji A72_i)(t3 At^—I Ayj). 

A direct calculation yields 

a + f3 


Bn = To - N - 1 - 


{a'^-/3^){2N + 2 + a + l3) 


1 


1 


and 


Un = 


a + 13 + in+ 2 a + f3 + 2n 
4,n{n -\- a){n -|- j3){n -I- a -I- (3){n -I- fV + 1 -I- a -I- j3){N -|- 1 — n) 


(2n + a + fi — l)(2n -I- a -I- /3)^(2n -|- a -|- ,5 -I- 1) 

It is seen that these recurrence coefficients correspond to the generic Hahn polyno¬ 
mials Qn{x',a,(3-,N); see [13] for the definition. 

The operators L and M defined as in dm and dSD) also generate a quadratic 
algebra. Let Z be the operator defined by 

Z = [L,M] = 2a;(l - x){2x - l)dl - 2Vi{x)d^ - Vo{x), 

where 

Vi (x) = (ts + 2a + 6 + 2j3)x^ — (ts -I- 6 -I- 3a -I- j5)x -I- a -I- 1, 

Vo{x) = {a + 13 + i + To){{a -\- (3 + 2)x - a - 1). 

It is verified that L, M and Z satisfy the commutation relations 

[Af, Z\ = 2 M'^ — 2(t3 -|- 2to)M — A L 2to(to + T 3 ) -|- 2(a -I- l)(/3 -b 1), 

[Z, L] = 2{L, M} - 2{t3 + 2to)L - (a + /3)(a + /3 + 2)M 

-b (a -b /3)(To(a -b /3 -b 2 ) -b T 3 (a -b 1 )). 

These relations correspond to the quadratic Hahn algebra |S], which is a special 
case of the Racah-Wilson algebra dm- 

7. Conclusion 

In this paper, we have explained from an algebraic standpoint how the tridiago- 
nalization of the hypergeometric operator leads to the Wilson polynomials; that is 
we have shown how the operators involved generate the Racah-Wilson algebra. 

It would be of interest to enlarge the scope of our study to cover g-orthogonal 
polynomials and their q —1 limits. In this context, it would be desirable to arrive 
at the Bannai-Ito algebra by the tridiagonalization approach and to determine 
how the Bannai-Ito polynomials can be obtained from lower families of orthogonal 
polynomials of the —1 hierarchy. 
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